Phase Diagram of Na^CoC^ Studied by Gutzwiller Density Functional Theory 
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The ground states of Na^CoC^ (0.0 < x < 1.0) is studied by the LDA+Gutzwiller approach, 
where charge transfer and orbital fluctuations are all self-consistently treated ab-initio. In contrast 
to previous studies, which are parameter-dependent, we characterized the phase diagram as: (1) 
Stoner magnetic metal for x > 0.6 due to a\ g van-Hove singularity near band top; (2) correlated 
non-magnetic metal without e' g pockets for 0.3 < x < 0.6; (3) e' g pockets appear for x < 0.3, and 
additional magnetic instability involves. Experimental quasi-particle properties is well explained, 
and the a\ g -e' g anti-crossing is attributed to spin-orbital coupling. 

PACS numbers: 71.20.-b, 71.27. +a 



Transition-metal oxides have complex phase diagrams 
due to the interplay between the charge, spin and orbital 
degrees of freedom. Among them Na^CoG^ is a typi- 
cal system showing doping-dependent phase control [l(. 
It has been found experimentally as non-magnetic (NM) 
metal for Napoor side, while Curie- Weiss metal for Na 
rich side 0, 0], and A-type (layered) anti-ferromagnetic 
(AF) state for x ~ 0.75 [3J. In addition, superconduc- 
tivity is discovered for hydrated NazCoGvyH^O (x ~ 
0.35) and charge-spin-orbital ordered states are sug- 
gested for x = 0.5 due to Na ordering The rich prop- 
erties of Naa;Co02 attract much of the research interests 
due to not only its potential applications, but also the 
challenging theoretical issues in this system generated by 
both the multi-orbital nature and strong e-e correlation. 

Naa;Co02 is crystallized in planar triangle lattice, with 
each Co site being coordinated by edge-shared oxygen- 
octahedron. The e g states are about 2eV higher than t^g, 
and the Fermi level is located within the Co-i2g multi- 
plet, which splits again into one a\ g and two e' g orbitals 
under trigonal crystal field. For Na concentration x, the 
effective number of t2 g electrons per Co is given as 5+x, 
and thus the low energy physics here is dominated by the 
multiple orbits (ai g +e' g ), where charge, spin and orbital 
degrees of freedom are all active. The rigorous compu- 
tational tools for such systems are still lacking, and the 
observed rich phenomena remains far from even qualita- 
tive being understood. The main controversial issues are: 
(1) for x = 0.3, are there any Fermi surface pockets for 
the e' g band? These pockets are predicted by LDA (local 
density approximation) calculations [f| but not observed 
by ARPES 0. (2) Is the x > 0.5 side more "correlated" 
than x < 0.5 side, as suggested by the Curie- Weiss be- 
havior for the Na rich side? It is expected from simple 
band picture that x = 1.0 end compound is a band insu- 
lator rather than Mott insulator. 

Both LDA and LDA+U methods fail for such system 
due to the insufficient treatment of electron correlation. 
This is why many issues look controversial following the 
LDA pictures. For instances, LDA predicts ferromag- 



netic metal as the ground state for the whole doping 
region, and LDA+?7 even enhances the tendency to be 
ferromagnetic [8]; the band width obtained by LDA or 
LDA+[/ is about two times larger than what observed by 
ARPES 9]; the e' g pockets problem as mentioned above; 
and etc. To treat the electron correlation more precisely, 
the Gutzwiller [1(3] and DMFT (dynamic mean field the- 
ory) [ll| approaches has been adopted, where fluctua- 
tion effects are included. However, those studies are only 
focused on the x — 0.3 compound using tight-binding 
Hamiltonian extracted from LDA, and conflicting results 
are drawn due to different parameters [T3 |. 

In this paper, we show that the above mentioned the- 
oretical challenging issues of this multi-orbital correlated 
electron systems can be well studied by using the re- 
cently developed LDA+Gutzwiller method [13j, which 
keeps the parameter-free character of density functional 
theory (DFT) and includes all possible charge transfer 
and orbital fluctuation effects self-consistently. As the 
results, a phase diagram for the whole doping region is 
constructed, which removes away most of the above men- 
tioned controversial issues, and suggests that physics here 
dominates by the doping-dependent orbital, charge and 
spin fluctuations. 

To overcome the problem of LDA, the common proce- 
dure of LDA+U and DMFT schemes [HI are to draw out 
from the LDA Hamiltonian the interaction terms for the 
localized orbitals, such as the 3d or 4/ states, and then 
treat the interaction Hamiltonian explicitly in a proper 
way (beyond LDA). The total Hamiltonian reads: 
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where \ia) a are a set of local orbitals with spin index a 



and occupation number n ia for lattice site i; the U, U' 
and J gives the intra-, inter-orbital repulsive interaction 
and Hund's exchange coupling, respectively. The H^ c is 
the double counting term from LDA, where interaction 
strength U and J are averaged over orbitals. 

If the interaction term is treated by the Hartree-like 
scheme (LDA+t/), the correction over LDA is a set of en- 
ergy shift of the local orbitals, leaving the kinetic part un- 
changed. This is fine if the fluctuation effect is not strong, 
but will fail in opposite case, such as the Na^OoC^ sys- 
tem studied here. In this sense, the DMFT method, in 
which frequency-dependent self-energy is properly com- 
puted, is much better than LDA+U. However, due to the 
heavy computational cost for multi-orbtial systems, the 
current DMFT studies " 
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[12| are all applied to tight- 
binding Hamiltonians extracted from LDA without full 
charge density self-consistency. This is insufficient if the 
charge, spin and orbtial degrees of freedom are all active 
as discussed above for Na^CoC^. 

In the LDA+Gutzwiller approach [13], the Gutzwiller 
wave function \^g) = P\^o) {P is a projection to many- 
body configuration) is used instead of single slater deter- 
minate wave function \^o). The orbital, charge and spin 
fluctuations can be included by the multi-configuration 
nature of the Gutzwiller wave function. As the results, 
a set of orbital-dependent kinetic energy renormaliza- 
tion factor Z a are obtained for the correlated states in 
addition to the on-site energy shift. Unlike the previ- 



ous Gutzwiler or DMFT studies [13, EH, 0, here all 
charge transfer processes, crystal field and orbital fluctu- 
ations are self-consistently treated within the framework 
of DFT, which allows for the accurate computation of 
ground state total energy. 

We use the plane-wave pseudo-potential method, and 
choose the Co-3<i wannier functions as the correlated lo- 
cal orbitals. The atomic-limit convention U = U' + 2J 
is followed. How to determine the value of U and J is a 
common problem for LDA+J7, DMFT and present meth- 
ods, and no unified way is established yet. Nevertheless, 
reasonable estimations have been done for U = 3.0 ~ 
5.0eV and J ~ l.OeV for Naa;Co02 system following the 
literatures 13, U|. Instead of using single fixed U, var- 
ious values have been studied, and our qualitative re- 
sults are not changed, as shown below. In addition, since 
our main purpose is to establish a general picture for 
the physics of Na^Co02, the structure differences among 
different doping x are neglected, and the Na doping is 
treated by virtual crystal approximation (therefore the 
charge-ordered states with Na ordering at particular dop- 
ing, which is interesting but not our purpose here, is out 
of the phase diagram). 

Fig.l shows the phase diagram computed for the whole 
doping range 0.0<a;<1.0. The solid lines and the dashed 
lines represent the stabilization energies of FM state and 
layer-type AF state relative to NM solution, respectively. 
The intra-plane AF state is hard to be stabilized due 
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FIG. 1: The calculated stabilization energies of FM (solid 
lines) and layer-type AF (dashed lines) states with respect to 
NM state for Na^CoC^ over the whole doping range 0.0 < 
x < 1.0. (a) Results obtained by LDA; (b) LDA+Gutzwiller 
results obtained using U = 3.0 eV and J = 1.0 eV; (c) 
LDA+Gutzwiller results obtained using U — 5.0 eV and 
J = 1.0 eV. 



to geometrical fluctuation of triangle lattice. The LDA 
(shown in Fig. 1(a)) gives magnetic ground states for all 
doping x, which are inconsistent with experiments. In 
contrast, the phase diagram by LDA+Gutzwiller has 
three distinct regions, which can be understood as the 
consequence of competition among crystal field splitting, 
inter-orbtial charge and spin fluctuation, as discussed in 
the following parts. As shown in Fig. 1(b) and (c), the 
features of phase diagram are qualitatively the same us- 
ing [/=3.0 eV or [7=5.0 eV. 

Phase II: correlated non-magnetic metal (0.3<x<0.6). 

First of all, the NM state is now correctly predicted 
for this region. It is known that LDA overestimates the 
tend ency to be FM for several systems, such as ruthen- 
ates [16|, and this artifact is even enhanced by LDA+C7. 
The physical reason is that correlation effect is not prop- 
erly treated in LDA and LDA+C7, but it is well included 
in the present formalism. In Fig. 2 we summarize the 
properties of the NM solutions for the whole doping 
range. For x larger than 0.3, as shown in Fig. 2 (a), the 
e' g bands are fully occupied and the inter-orbital fluctu- 
ation (defined as F and S, see caption of Fig. 2) is weak, 
indicating a effective single band system. However with 
x approaching the phase boundary around x c = 0.3, a 
crossover to multi-band behavior has been detected from 
the strong inter-orbital spin and charge fluctuation as 
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FIG. 2: Results for NM state as the functions of dop- 
ing x with U=5. OeV and J=1.0eV. (a) Occupation num- 
bers; (b) Band renormalization Z 2 -factors, inter-orbit charge 
fluctuation F=(h aig n e ' ) — {nai g ){n e ' }, inter-orbit spin fluc- 
tuation S={Sa lg S*,) _ (&i g )(8t'h and (°) Level splitting 

A = £„,„ ~ EJ . 



shown in Fig. 2(b). Such fluctuations in low doping area 
is induced by the Hund's rule coupling, which favors even 
distribution of electrons among different orbitals with 
same spin. Due to the presence of the correlation effect, 
the quasi-particle band width (kinetic energy) is renor- 
malized by factor Z\ which is about 0.5 (0.7) for the 
a\ g (e' g ) state at £=0.3 (as shown in Fig. 2). The same 
amplitude of renormalization is reported by ARPES 0] • 

Secondly, it has long been a controversial issue whether 
the e' g states cross the Fermi level? To answer this ques- 
tion, the correlation renormalized level shift is crucial 
(as shown in Fig. 2(c) the a\ g -e' g level splitting A is much 
renormalized compared to LDA results). Unfortunately 
after including the correlation effect, two studies have 
been done, and conflicting results are drawn 13, HH for 
£=0.3. It was recently pointed out by Marianetti et al. 
that the controversial is due to the different choice of 
crystal- field splitting in their tight-binding model [12 |. 
To go further, in our studies, not only the crystal-field 
is treated parameter-free, but also the full charge self- 
consistency is achieved. To see the difference, we per- 
formed one-loop calculations (i.e. the charge-density is 
fixed to the LDA value and only Gutzwiller wave func- 
tions are optimized), then Marianetti's results are recov- 
ered, i.e, the e' g pockets are not present for £=0.3. In 
addition, we also found that even for a;=0.2 the e' g pock- 
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FIG. 3: The calculated band structure of Naa;Co02 for 
(a)a;=0.35; (b)x=0.5; (c)x=0.75 in the NM state, (d) is the 
projected density of states (PDOS) for ir=0.8 in the A-type 
AF state. U=5.0 eV and J=1.0 eV are used, and the spin- 
orbital coupling is included. 



ets are not present in this one-loop calculation. How- 
ever, after including the charge density self-consistency, 
the renormalization of level splitting are suppressed (see 
Fig. 2(c)). As the results, we found that a;=0.3 is the 
critical point, i.e. e' g pockets are absent for x >0.3 but 
present for x <0.3. Our calculated quasi-particle bands 
(as shown in Fig. 3 and discussed below) can be well com- 
pared with ARPES. Since with the interactions we used 
(£7=5. OeV and J=1.0eV) the static Hartree Fock shift is 
zero , the renormalized energy level shift here is fully 
contributed by the fluctuation effect. From Fig. 2(c), we 
found that the A is peaked in the crossover region indi- 
cating that the inter-orbital fluctuation is the main rea- 
son for the renormalization effect on the energy levels. 

Phase III: weakly correlated Stoner metal for x > 
0.6. Clearly seen from the calculated density of states 
(Fig. 3(d)), a sharp van-Hove singularity (VHS) is present 
near the a\ g band top-edge (due to the flat dispersion). 
For the Na rich side, the Fermi level is shifted close to 
the VHS, the Stoner instability make the system FM 
(in-plane). This conclusion is supported by the follow- 
ing facts for this region: (I) LDA works qualitatively 
well; (2) the magnetic solution only weakly depends on 
interaction strength U and J (see Fig.l); (3) calculated 
spin and charge fluctuation are all weak (Fig. 2(b)). It 
is therefore suggested that strong correlation is not the 
driving force for the magnetic state, instead the VHS 
is responsible. Furthermore, we correctly predict that 
the A-type AF state is more stable than FM state (see 
the difference between solid and dashed lines in Fig.l). 
From Fig. 3(d) for £=0.8, it is seen that spin moment 
mostly comes from the a\ g states, which aligns towards 
the z direction (inter-layer) of the crystal. According to 
the Goodenough-Kanamori rule, the exchange-coupling 
along z direction is dominates by AF super-exchange, 
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which stabilizes the A-type AF state. In fact, the esti- 
mated inter-layer exchange coupling (from the total en- 
ergy difference between FM and A-AF solutions) is about 
J c =3.0 meV, in good agreement with experimental re- 
sults Q . From itinerant FM theory, the size of spin mo- 
ment may change with raising temperature, but in the 
presence of strong VHS near band edge, this possibil- 
ity is prevented by the sharp density barrier. Indeed we 
found that the calculated moments of FM and A-type AF 
solutions are the same. Therefore, the Oi ff moment will 
behavior like localized spin as observed experimentally. 
In other words, the Curie- Weiss behavior in this region 
does not necessarily suggest the enhanced correlation. 

Phase I: magnetic correlated metal: The main differ- 
ence between phase I and II is that the e' g band start 
to go cross the Fermi level and the e' g hole pockets are 
present. Therefore phase I is effectively a multi-band sys- 
tem. In the meanwhile, the magnetic instability is recov- 
ered, and the system is stabilized in FM ground state. 
Several points should be addressed here to understand 
this phase: (1) in contrast to phase III, the FM state is 
slightly more stable than A-type AF state; (2) both e' g 
and a\ g contribute to the spin moment, in other words, 
the stabilization of FM state is due to the enhancement 
of inter-orbit (rather than intra-orbit) spin fluctuation 
(Fig. 2(b)), which is induced by the Hund's rule coupling; 
(3) strong correlation (rather than VHS) is responsible 
for the magnetic instability. It is interesting to note that 
the critical point x ~ 0.3 (boundary between phase I 
and II) is close to the doping level where superconduc- 
tivity was observed. The experimental information for 
this region is quite limited and not conclusive due to the 
difficulty of sample preparation, our prediction should be 
evaluated by future experiments. 

Finally, we show systematically in the Fig.3(a)-(c) the 
calculated quasi-particle band dispersion for x=0.35, 0.5, 
and 0.7 after including the spin-orbital-coupling (SOC) 
effect. The overall picture can be nicely compared to 
ARPES data 0- In particular, (1) the band width renor- 
malization around factor of 2 is now obtained; (2) the 
e' g Fermi surface pocket is absent; (3) the a\ g -e' g anti- 
crossing along r — K line is nothing but a effect of SOC, 
and the gap around 0.1 eV is comparable to experimental 
data M- 

In summary, using the recently developed 
LDA+Gutzwiller method, we are now able to cal- 
culate the ground state total energy of correlated 



multi-orbital systems from ab-initio after taking into 
account the orbital fluctuation. The calculated phase 
diagram of Na x Co02 establishes a general understand- 
ing for the physics behind. Most of the discrepancies 
between experiments and previous theories, such as the 
e' g pocket, a\ g -e' g anti-crossing, Curie-weiss behavior for 
x >0.5, are self-consistently understood. Three distinct 
phase regions are identified, which is instructive to 
future experiments. 
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